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^ ' Abstract 

X> . 

In general, while obtaining the probability density function of sums and products of shifted 
random variables, ordinary analytical methods such as Fourier and Mellin transforms tend 
to provide integrals which cannot be expressed in terms of ordinary Meijer G and H func- 
tions. This way, the need of defining new functions which easily enable one to write such 
integrals in a closed-form is inherent to the development of this area of statistical sciences. 
By generalizing the Mellin transform which defines the H function, a new function is estab- 
lished. A direct application of the so-called I is discussed while developing the probability 
i -^h , density function of the sum and the product of shifted generalized gamma random vari- 

ables. Important special cases of the I and their applications in science are also discussed 
in order to show the applicability of the function hereby defined. 

Keywords: generalized gamma distribution, shifted random variables, Mellin transform, 
Fourier Transform, generalized H function. 
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The use of shifted distributions in the modelling of practical situations has grown considerably 
over the last years. Such growth is deeply related to the possibilities that arise by inserting 
a shifting - or translational - parameter into the definition of a given random variable. For 
example, the possibility of positioning the center of mass of a given distribution along its 
support enables a better description of the data analyzed. Also, while modelling measured 
data, shifting parameters make it possible to correct some experimental limitations such as 
measurements above/below a given threshold. 

Even though shifted distributions have become more frequent in the scientific community, 
standard analytical functions and methods which would be capable of dealing with this kind of 
distributions fail in their intent. This way, a new approach to this problem is required. 

In the present paper, a new function, hereby named I function, which generalizes the well- 
known H function is defined. 

In short, this new function changes the usual gamma functions in the contour integral 
representation by H functions themselves. This way, the / function contemplates as special 
cases the H function the I function [5], the Y function [8J, the Generalized Upper Incomplete 
Fox H function |7j and other hypergeometric-type special functions. 

In order to show the applicability of the I function, the latter is used to represent the 
distributions of the sum and the product of generalized shifted gamma random variables. This 
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result could not be achieved by means of the existent functions, which ensures the need of the 
new function defined in the present paper. 

In the next section, the definition of the H function and its Mellin transform are presented 
which are used in the definition of the / function. 



2 The H- function 



The H - function (see |6J,[4J and [3j) is defined as an contour complex integral which contain 
gamma functions in their integrands by 



Tjm.n 
H p,q 



(ai,Ai), (a n ,A n ), (a n+1 ,A n+1 ), (a p ,A p ) 

(bi,Bi), (b m ,B m ), (b m+ i, B m+ i), (b q ,B q ) 



2tu 



Yl r{bj + B jS ) JJr(l - a 3 - A 3 s) 
j=i 5=i 

q p 

J] T(l-b 3 -B 3 s) J] T(a j + A jS ) 

j=m+l j=n+l 



-z s ds. 



(1) 



where Aj and B 3 are assumed to be positive quantities and all the a 3 and b 3 may be complex. 
The contour L runs from c — ioo to c + ioo such that the poles of Y(b 3 + B 3 s), j = 1, . . . , m lie 
to the left of L and the poles of T(l — aj — Ajs), j = 1, . . . , n lie to the right of L. 
The Mellin transform of the H -fuction is 
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Given these, one shall proceed to define the generalized H function. 



(2) 



] | T(l - b 3 - B 3 s) Y[ T{a 3 + Aj S ) 

j=m+l j=n+l 



3 The Generalized H function: The / function 

The generalized H function, hereby named I function, can be defined as a contour complex 
integral which contain H functions in their integrands. For simplicity and in order to address 
a few problems of interest, the following simplified definition can be given 
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in which T(s) is the Mellin transform of the new function and can be explicitly given as: 
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where Aj^ and Bj^ are assumed to be positive real quantities, the a^fc, dj t k, bj,k, bj,ki II.,-, 
Wj, Tj, 7., p.,-, CTj, a,, Aj, Qj, G, ?7j, j = 1, ...,777, are real numbers. The contour L runs 
from c — ioo to c + zoo, where c is a real number, and exists in accordance to Mellin inversion 
theorem, taking into account all the singularities. 

Based on the definitions above, direct applications of the new function hereby introduced 

m 2 1 

are shown in the next section. The new function can be further generalized by changing H 3 ' 2 
to H™^ 1 . This latter consideration is discussed in Section 6 of the present paper. 

4 Generalized Gamma Distribution and Distribution of its 
Product and Sum 

In the present section, a direct application of the generalized H function is developed. Let the 
probability density function of the generalized gamma distribution be given as: 

OS 

/ (ar) = JP-L {x - ^f- 1 e-fa-nr, x > ^ «, /3, 7 > (5) 
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One may notice that, in (J3]), by setting a = 7 = k; fi = and = X~ k , the Weibull 
distribution with shape parameter k and scale parameter A is retrieved. On the other hand, in 
(jSJ), when a — 7 = 1; \i — and — A, the Exponential distribution with rate parameter k is 
obtained. Also, by setting /i = 0,a = 7 = 2 and = 1 /2a 2 , the Rayleigh distribution with 
parameter o is recovered. By taking a = 3, /x = 0, 7 = 2 and = l/2a 2 a Maxwell-Boltzmann 
distribution with parameter a is obtained. Finally, by setting /x = 0, a = k, 7 = 1 and = A, 
the generalized gamma distribution reduces to an Erlang distribution with shape parameter 
k and rate parameter A. It is worth noticing that other gamma-type distributions may be 
expressed as special cases of (J3J). The Generalized H function can be used to express (j3j) as 
follows: 

Theorem 1. The probability density function of a generalized gamma random variable is given 
as: 
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Proof. Consider the Mellin transform of the probability density function of a generalized gamma 
distribution. By means of (j^J), one may get: 
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On the other hand, consider the Mellin - Barnes representation of the exponential function 
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This way, by inserting (jSJ) into ([7]) , the Mellin transform given in (J7J) can be written in terms 
of the H function as: 
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M[f(x)}(s)= 1 ^fi a - 1+s H 2 3 j 

This way, by means of the inverse Mellin transform theorem, the representation follows 
from Q. □ 

N 

Theorem 2. Consider the random variable X = Y\Xi in which X{, i = 1, .., N are independent 

generalized gamma random variables. This way, the probability density function of X is given 
by: 
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Proof. It is known that the Mellin transform of the distribution of the product of independent 
random variables is the product of the Mellin transforms of each variable [6J, this way, the 
Mellin transform of the distribution of the product of N independent generalized gamma random 
variables is easily given by means of fl9]) as: 

N ^ 

M[fx(x)}(s) = l[ 1 ^^- l+s Hli 

i=i 

in which the subindex j indicates the parameters of each distribution. By means of ([3]) and 
01]), the representation in (ITO]) easily follows. □ 

It is also worth noticing that the quotient of generalized gamma random variables can be 
easily obtained by making the substitution s = 2 — s in Q) for the random variable which is in 
the denominator of the ratio [6J. Due to the simplicity of the procedure, its full development is 
not presented in the present paper. 

N 

Theorem 3. Consider the random variable Y = J^Xj in which X i} i — 1, .., N are independent 

3=1 

generalized gamma random variables. This way, the probability density function of Y is given 
by: 
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Proof. In order to obtain the probability density function of the sum of N independent shifted 
generalized gamma random variables, consider the Laplace transform of the probability density 
function given in (jSj) 



L[f ■(*)](*) = JT z 
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On the other hand, consider the Mellin - Barnes representation of the exponential function 
e -/3y 7 ? given in (jSJ). 

This way, by inserting dHJ) into (ITS"]) , the Laplace transform given in (ITS"]) can be written in 
terms of the H function as: 
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It is known that the Laplace transform of the distribution of the sum of independent random 
variables is the product of the Laplace transforms of each random variable [6], this way, the 
Laplace transform of the distribution of the sum of N independent generalized gamma random 
variables is easily given as: 
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By means of fl3]) and (j3J), the representation in (TT2l) easily follows. 



(15) 
□ 



The formulas obtained for the sum of N independent generalized gamma random variables 
easily enable one to obtain the distribution of the linear combination of random variables of 
this type. This can be achieved by noticing that, if X is a generalized gamma random variable 
with parameters a, (3, 7 and /x, for a given constant value A, the random variable AX is also 
a generalized gamma random variable with parameters a, jS/A 1 , 7 and Afi, respectively. 

One may also notice that the results presented in the present section are a generalization 
of the ones in [7] and [8J. Also the new results are entirely dependent on the definition of the 
I function. 



5 Special Cases of the / function 

In the present section, a few special cases of the generalized H function are discussed. 



5.1 Standard H function 

The key equation to relate the new function hereby defined to the standard H function is the 
contour integral representation of the gamma function, given as: 
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On the other hand, by setting aj = £ • = fjj = Tj = IL, = Aj = p ,. = 0, 7^ = Oj = = 7Tj = 
= 1, and taking into account ( ITT)) the first produtory present in the numerator of (120)) is 

recovered. The other produtories are easily obtained by similar procedures, ultimately showing 

that the H function is a special case of the I function. 

5.2 I - function 

In [5J, a generalized hypergeometric function has been defined. The main advantage of the so- 
called I function is that powers of gamma functions are considered inside the Mellin transform 
used to define such a function. Explicitly, in [5] it has been defined the I function as the inverse 
Mellin transform of the following function I(s): 
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This way, by setting aj = (j = rjj = Tj = EL, = Aj = pj = 0,jj = Qj = f3j = Wj = 
l,(fj = Bj, and taking into account (TTT)) the first produtory present in the numerator of (118)) is 
recovered. The other produtories are easily obtained by similar procedures, ultimately showing 
that the I function is a special case of the I function. 

It is interesting to notice that recently the I function has found important application 
in the study of wireless communication systems. In [2] it has been derived the probability 
density function and cumulative distribution function of the sum of L independent but not 
necessarily identically distributed Gamma variates. Their result is applicable to the modelling 
of the output statistics of maximal ratio combining (MRC) receiver operating over Nakagami- 
m fading channels. The results are derived in terms of Meijer G function, H function and I 
function. Following a similar approach as in [2] , in [1] it has been derived the probability density 
function and cumulative distribution function of the sum of L independent but not necessarily 
identically distributed squared rj — p, variates. As in the case of [2J, the results are applicable to 
the output statistics of maximal ratio combining (MRC) receiver operating over rj — p fading 
channels. The I function has a major role in the development of the results of both [2] and [TJ. 



5.3 Y - function 

In a recent paper [8], it has been proposed a generalization of the H function. In short, the so- 
called Y function generalizes the Mellin transform of the H function by changing the ordinary 
gamma function in ($2§ by Tricomi hypergeometric functions, defined as the following contour 
integral: 

U (a, 6, z) = — —— - — - — [ T(s)r(l - b + s)r(a - s)z- s ds, (19) 



r (a) T (a - b + 1) 2ni J c 

This way, by pre-multiplying the Tricomi hypergeomtric functions by a constant raised to 
a linear function of s, in [8] it has been defined the Y function as the inverse Mellin transform 
of the following function Y(s): 
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One may notice that (120]) is misprinted in [8j. 

In order to present the relation between both functions, by means of ([I]), it is easy to notice 
that: 



U (a, b, z) = H 2 A 
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Finally, by setting aj = (j = rjj = Fj = Uj = pj = 0,/3j = 7^ = Bj,Aj = f3j,Qj = 
<pj + bj — = Wj = 1 and by means of EU the first produtory present in the numerator of 
f l20|) is recovered. The other produtories are easily obtained by similar procedures, ultimately 
showing that the Y - function is a special case of the I function. 

The Y function has found interesting applications in wireless communication systems. In 
special, |S] obtained the PDF and CDF of a shifted gamma random variable in terms of the 
Y function. By means of their results, it was possible to simulate the outage capacity of a 
multicarrier transmission system through a slow Nakagami-m fading channel. 

5.4 Generalized Upper Incomplete Fox H function 

Also in recent paper [7j, it has been proposed another function which generalize the H function. 
The so-called Generalized Upper Incomplete Fox H function generalizes the Mellin transform of 
the H function by changing the ordinary gamma functions in (j2J) by upper incomplete gamma 
functions, which may be defined by means of the following contour integral: 

1 ' 2m] c r(l + s) 1 ' 

In [7J it has been defined the Generalized Upper Incomplete H function as the inverse Mellin 
transform of the following function UI(s): 
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In order to present the relation between both functions, by means of ([I]), it is easy to notice 
that: 



T (a, z) = H 2 3 j 



(0,0), (1,1), (1,0) 

(o,i),(M) 



(24) 



Thus, by setting aj = (j = r/j = Fj = Uj = Aj = pj = 0, 7^ = Bj, Qj = f3j = tcj = aj = 1, 
and taking into account f )24p the first produtory present in the numerator of ( I23p is recovered. 
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The other produtories are easily obtained by similar procedures, ultimately showing that the 
Generalized Upper Incomplete Fox H function is a special case of the / function. 

As in the cases of both I and Y function, the Generalized Upper Incomplete H function has 
found useful applications in the modelling of wireless communication systems. In [7J the PDF 
and CDF of the product of shifted exponential random variables has been obtained in terms of 
the Generalized Upper Incomplete H function. By applying such results, it has been possible 
to model the outage capacity of a multicarrier transmission system through a slow Rayleigh 
fading channel. 



6 Further generalization of / function 

A generalization of the / function, the J function, can be defined as follows: 
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in which T(s) is the Mellin transform of the new function and can be explicitly given as: 
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where, as discussed in Section 3, Ajj. and Bj ^ are assumed to be positive real quantities, 
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1 m are real numbers. 



The contour L runs from c — ioc to c + ioo, where c is a real number, and exists in accordance 
to Mellin inversion theorem, taking into account all the singularities. The parameters mj,rij,Pj 
and <L, j = 1, ...,m, are positive integers following H function's definition. It is important 

to notice that the braces sign used in the definition of the I function denotes a vector whose 
components are to be written as the generic case shown. 



7 Conclusions 

In spite of the growth in the use of shifted random variables, analytical functions and methods 
are unable to give closed-form representations for sums and products of this type of random 
variables. 

In the present paper, in order to approach this issue, a new / function, is defined. This new 
function generalizes the H function by changing the gamma functions present in its contour 
integral definition by H functions themselves. This way, not only the H function but also the 
I, Y and Generalized Upper Incomplete Fox H functions are shown to be special cases of the / 
function. 

The new function has been successfully applied to obtain the distributions of the sum and the 
product of generalized shifted gamma random variables, showing its value both in theoretical 
and practical backgrounds. 
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